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Abstract. We present some results of an ongoing study of the long-term 
spectroscopic variations of AG Dra, a prototypical eruptive symbiotic system. 
We discuss the effects of the environment and orbital modulation in this system 
and some of the physical processes revealed by a comparison with the nova 
outburst of the symbiotic-like recurrent nova V407 Cyg 2010. 

Key words: Symbiotic stars; stars- individual (AG Dra, V407 Cyg); stars: 

novae; physical processes 

1. INTRODUCTION 

The symbiotic binary AG Dra certainly needs no introduction to readers of 
these proceedings and we refer to Leedjarv’s report in these proceedings for the 
relevant background and the excellent introductions in Leedjarv et al. (2004), 
Viotti et al. (2007), Galis et al. (2007), and Gonzalez- Riestra et al. (2008). Pho- 
tometric monitoring of this system extends now over nearly a century with the 
highest precision spanning more than 40 years from both amateurs and profes- 
sional astronomers (see Munari, these proceedings, and Skopal et al. (2009)). We 
concentrate, instead, on the line profiles and flux measurements. We have used 
new and archival optical grating, grism, and echelle spectra from Asiago (Pennar 
and Ekar), Loiano, the Telescopio Nazionale Galileo (TNG), Ondrejov, and the 
Nordic Optical Telescope (NOT). Part of these data have already been published 
(Shore et al. 2010) and we refer the reader to that paper for details covering the 
2006-2009 part of the observations. The full details of the data set are available 
in Genovali (2010). Here we outline some of the results from our study of AG Dra 
that began with the major outburst of 2006-2008 and has since expanded to cover 
the last 30 years (Genovali et al. in preparation). 
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Fig. 1 . The Balmer line variations during the 2006-2008 outburst of AG Dra, corrected 
for variations in the continuum flux as described in Shore et al. (2010). Square: Ha; 
cross, H/3; asterisk: H 7 ; diamond, H<5. The outburst covered the orbital interval 5. p 88- 
7. p 15 with the peak fortuitously occurring near superior conjunction of the giant. 

2. EMISSION LINE FLUX VARIATIONS 

We begin with two summary plots. Figure 1 is the variation of the Ha through 
Hd line equivalent widths, corrected for the continuum variations, from the major 
outburst of 2006 of AG Dra. Figure 2 shows the relative He I, He II 4686A, and O 
VI 6825A variability during the 2006 outburst. We draw the reader’s attention to 
the bifurcation in the last panel that may distinguish the hot and cool outbursts; 
those from the start of the 2006 eruption are the upper branch for He I. As noted 
by Gonzalez-Riestra et al. (2008) and Skopal et al. (2009), there is a saturation 
in the He II emission level during the major eruptions. This could be due to the 
increased optical depth of the wind of the hot component. Viotti et al. (2007) have 
already conjectured that a burst may have two stages depending on its intensity: 
a hot event may, if it continues long enough, become a cool outburst. 

Passing to the long haul, we show in Fig. 3 the variations in the equivalent 
widths of the main emission features and the UBV photometry for about two 
decades. 1 It is an impressively rich phenomenology (see also Galis et al. 2007). 
Informed by the dense coverage of the 2006 outburst, going backwards in time to 
the major outburst of 1995, it appears that the pattern of Raman 6825A feature 
variation was the same. A steep rise from what was likely the disappearance of the 
line. A relatively minor outburst, about three years later, shows the same pattern. 
Leedjarv reported at the meeting that the same thing appears to be happening 
again with the most recent (2010) minor outburst. Since the variations of this 
line are directly linked to the down-conversion of O VI 1032A photons, there are 
two possible, likely linked inferences: the column density in the wind around the 
white dwarf (WD) could have changed and/or the intensity of the resonance lines 
could have changed. The first implies a variability of the mass loss rate from the 
red giant (RG) or a change in the ambient density around the WD. The second, 
instead, is something intrinsic to the WD mass loss, a change in the optical depth 
of the O VI lines implying a variation in the wind from the degenerate component. 
Two close peaks around day 1000 (Fig. 3) show the same Raman line variation 


1 The displayed data are from Leedjarv et al. (2004) (star), Tomova & Tomov (1999) (pen- 
tagon) and spectra from Loiano (plus), Asiago/Ekar (asterisk), Catania (diamond), TNG (cross), 
Asiago/Pennar (filled circle). In the lower panel: U (filled circle), B (plus), and V (circle) from 
Skopal and collaborators and the AAVSO. 
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Fig. 2. Variations of the He I, He II, and O VI Raman 6825A line during the 2006- 
2008 eruption from the Loiano and Ondrejov spectra. Top left: He I 6678A (square). He 
I 7065A (cross), and He II 4686A (diamond) variations in time. Top right, He I 6678A 
vs. He II 4686A. Bottom left, He II 4686A vs. O VI Raman 6825A; bottom right: He I 
6678A vs. O VI Raman 6825A. Note, in the last plot, the bifurcation. 

eventually reaching the same maximum. If these were so-called hot outbursts, the 
O VI line would have turned less optically thick than in the cool outbursts, and the 
feature would not have disappeared. If, instead, the WD mass loss was higher, and 
the covering factor complete, then the increase in the WD wind could have been 
so large that the O VI 1032, 1036A was smuthered with a consequent vanishing 
of the optical features (cf. Schmid (1995) and Harries & Howarth (1997) for a 
discussion) . 

The 2006 outburst confirmed what had, implicitly, already been reported in 
the literature, that the Raman lines can disappear (Munari et al. 2009). Having 
been seen on only one occasion, however, and lacking dense coverage it was not 
obvious that this is a generic feature of a certain magnitude of event. Instead, 
with 20-20 hindsight, it seems that the 6825 A feature flux curve observed during 
the last major outburst is actually typical, the slow recovery of the line relative 
to its sharp initial drop and its subsequent return to the quiescent state being 
recognizable in earlier events. The interval of disappearance was brief, about a 
week, and due to weather, season, and scheduling, other similar events might 
have been (barely) missed. The rise from minimum line strength is similar for 
all such large U outbursts. An important feature of the long term variations is 
the appearance of a sort of saturation level for the emission line strengths, in 
particular He II 4686A as a function of the U magnitude. From Leedjarev et al. 
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Fig. 3. The long term spectroscopic history of AG Dra. Equivalent width variations of 
(top to bottom) Ha, H f), O VI Raman 6825A, and He II 4686A emission lines, corrected 
for continuum variations, together with the UBV photometry. 
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(2004), it was already known that the emission varies in phase with the near UV 
but the relation was a simple power law up to U« 9 m . At brighter levels, the 
strict increase breaks down with a large dispersion, greater for He I 6678A than 
for He II 4686 A. The O VI Raman 6825 A feature actually peaks at U rjIO and 
then declines by about a factor of two for increasing brightness. The maximum 
equivalent width reached by the Raman 6825A line is -15A that is, however, rarely 
reached in the historical data. 

The optical He II lines, formed by recombination in the circum-WD region, 
are useful, but not unique, proxy measures of the state of the hot source. The 
recombination time is density dependent and can be of the same order as the 
interval between activity phases. It should be recalled that the environment is 
a mess given the many dynamical and radiative processes acting simultaneously. 
These lines convolve the ionization source variations with the response function 
of the environment, an unknown given the uncertainties regarding the mechanism 
producing the ultraviolet and X-ray variations. The Raman features are, instead, 
instantaneous measures of the FUV since they are optically thin and produced only 
by scattering. Using this, it is interesting that one recovers the same classification 
of hot and cool outbursts obtained from the direct measurement of the X-rays with 
ROSAT and XMM analyzed by Gonzalez-Riestra et al. (2008) and Skopal et al. 
(2009). Those events during which the Raman line fell below the quiescent floor 
and the U magnitude increased by the largest amount were judged cool, while 
the the two events that saw the line increase from the floor level and return to it 
after a short time were hot. The other intriguing feature is that during the hot 
outbursts the Raman line reaches a systematically greater flux. 

The H + region of the WD is always ionization bounded, we are not dealing 
here with an interstellar cloud, but it is dense and not without intrinsic line opac- 
ity that can affect the spectrum formation at other wavelengths. For example, 
changes in the Lyman series optical depths, that can be caused by variations in 
the FUV luminosity of the hot source or a change in the absorbing column den- 
sity, can only be extracted by detailed modeling of the line formation. Decreasing 
the density increases the ionized volume and decreases the Lyman series optical 
depth, changing the optical Balmer line intensities and profiles. Changing the 
mass loss rate from the hot star displaces the location of the stagnation point of 
the wind- wind shock and also may change the spectrum of hard photons available 
to further ionize the environment. This requires self-consistent, time dependent 
photoionization and hydrodynamic models of the sort described by Folini et al. 
(1998, 2003). 

Like the Raman events, there seems to be a characteristic time development 
of the U outbursts but these also indicate different natures. A change in the 
covering factor and/or the column density without an accompanying change in 
the luminosity of the hot source will produce only flux redistribution and the 
bolometric luminosity should remain constant. This is well known from the post- 
explosion stage of classical novae and is not far from correct for LBVs. In such 
cases, an increase in the mass loss rate from the hot star could have the same effect. 
Instead, as suggested by Sokoloski et al. (2006), if the outbursts are actually mini- 
TNR events (see Starrfield, these proceedings), then the bolometric luminosity 
should vary. At least the major event of 2006 appears to be the latter case. This 
is a further way to separate the sources of the variations. 
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Fig. 4. The equivalent width variations of the O VI Raman 6825A, and He I 6678A 
emission lines in AG Dra, corrected for continuum variations, phased on the orbital 
solution by Fekel et al. (2000). See text for details. 

3. ORBITAL MODULATION OF STRONG EMISSION LINES 

Line profile variations are direct probes of the environmental effects along the 
line of sight. To illustrate the complexity of the AG Dra phenomenology, we show 
some sample phased variations of the Ha and helium lines. Figure 4 shows the He 
I 6678A and O VI 6825A scaled flux variations plotted on the Fekel et al. orbital 
solution. No distinction is made in this plot between phases of activity from JD 
49030 to JD 54639 (Ondfejov and Loiano spectra only). For the Raman feature, 
the minimum line strength (< 0.3A) from the 2006 outburst happens to have 
occurred at orbital phase zero. This plot also suggests the presence of multiple 
periodicities in the line variations. While the Raman line shows little orbital 
modulation, there is a clear indication of such in the He I variations. Further 
data cuts, based on the activity state, show that the line fluxes from pre-selected 
quiescent stages are those with the strongest orbital modulation, as can be gleaned 
from Fig. 3. 

The Raman features are immune to self-absorption, that alters the permitted 
Balrner and helium profiles, and Rayleigh scattering that affects the UV lines. 
Electron scattering should not be important since it is local to the source, not to 
the conversion region. Instead, they should be sensitive to changes in the line of 
sight column densities of neutral hydrogen. 2 The modulation of the ratio of the 
two optical O VI Raman features may, consequently, be indicative of structural 
changes in the wind resulting from the hydrodynamic interaction of the winds 
from the two components (see Fig. 4 in Shore et al. 2010). The minimum in the 
line ratio happened in the interval between the two U maxima of the outburst and 
is, perhaps coincidentally, approximately at inferior conjunction of the WD. We 
urge observers to include, whenever possible, both features in future spectroscopic 
studies. In Figs. 5a, b and 6a, b we show sample high resolution line profiles at a 
few phases to illustrate the environmental effects. 

The He I line profiles show a nebular structure, the He II 4686, 5412A lines 
show, instead, a comparatively single peaked profile that remains at almost con- 
stant radial velocity, about -200 km s -1 , but whose symmetry changes. The vari- 
ations of the profile of the 1640A line, studied by Gonzalez-Riestra et al. (2008), 


2 A similar effect was discussed for the O I] 164lA line in symbiotics by Shore & Wahlgren 
(2010) as a probe of the neutral wind of the RG. In the absence of available UV measurements, 
the [O I] 6300A line, formed from the same transition cascade, may be an appropriate substitute. 
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Fig. 5a. Sample profiles for H a in AG Dra from the TNG spectra at several orbital 
phases. The dates are JD 51826, 53151, 53589, and 54731. The phases are according to 
the ephemeris of Fekel et al., (2000) with (j) = 0 P is RG inferior conjunction. 
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Fig. 5b. The Balmer H7 and H5 line profiles from the NOT spectrum on JD 53100, 
at orbital phase 0. p 23. 

are more extreme in displaying extended wings with maximum velocities of 2000 
km s -1 that do not appear to show orbital modulation. 
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Fig. 6a. Examples of the He I 5875A line profiles in AG Dra as a function of orbital 
phase from Catania, TNG, and Asiago spectra. The dispersions were 0.02A px -1 and 
0.15A px -1 . The time period is from JD 51826 to JD 55076 covering eight orbits. 
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Fig. 6b. Examples of the structure on the He I 6678A singlet line profiles; as in Fig. 
6a, from TNG and Asiago. 


Figure 7 shows a simple model of the line profile variations for Ho as a func- 
tion of orbital phase (see Dumrn et al. (1998) for an application to BX Mon). 
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Fig. 7. Model line profile variations with orbital phase. A synthetic Ha line, assumed 
to be a gaussian, was absorbed in the wind of the RG assuming a terminal velocity of 30 
km s -1 and a simple /3-law (see text). The mass ratio and orbital period approximately 
match BX Mon. Differential displacement (center-of-mass frame) produces the variability 
of the absorption feature. The emission line intensity was kept constant. 


The emission line from the ionized region around the hot star, assumed to be a 
gaussian for simplicity, is observed with orbital radial velocity shifts through the 
extended neutral wind of the RG that has a /3-type velocity law. In this case, a 
velocity amplitude of order 10 km s -1 for the RG and a mass ratio of 3:1 (RG:WD) 
were assumed to simulate BX Mon (see Dumm et al. (1998)). No quantitative 
information is intended here, just a warning that even in the visual part of the 
spectrum, environmental line absorption effects that are well known to affect ul- 
traviolet line profile studies can produce significant changes in the profiles. It then 
suffices to recall that in the real systems, the interactions of the winds of the two 
components, complex structures within the orbital plane, non-spherical outflows, 
and time dependent variations, combine to complicate the quantitative analyses of 
such variations. In AG Dra, for instance, the He II 1640A emission line intensity 
variations and profiles are clearly affected by both intrinsic changes and orbital 
modulation by the UV Fe-curtain and detailed models are needed to separate the 
various contributions. 

4. LESSONS LEARNED FROM THE 2010 NOVA ERUPTION OF V407 CYG 

One of the many things learned from the 2010 eruption of V407 Cyg, a symbi- 
otic that showed all the symptoms in this event of a symbiotic-like recurrent of the 
RS Oph variety (Munari et al. 2011, Shore et al. 2011a, b), is that the presence 
of broad wings on the emission lines, especially intrinsically asymmetric profiles, 
can also have hydrodynamic origins aside from steady outflows. The interactions 
in the environments of these systems are complex and very time dependent, espe- 
cially in the lower density parts of the circumstellar medium, as shown by two and 
three dimensional simulations. We show some examples of the late-time Balmer 
line profiles that allow a determination of the individual regions and their reaction 
to the expanding shock and photoionization changes caused by the post-eruption 
emission from the WD. The lines are asymmetric, a characteristic feature of an 
off-center shock propagating in the density gradient of the RG wind. Such profiles 
have, intriguingly, been observed on He II 1640A in AG Dra and other symbi- 
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Fig. 8. The Balmer line profile in the symbiotic-like recurrent nova V407 Cyg on 
2010 Jul. 15 (top) and 2011 Aug. 21 (bottom) (date of outburst 2010 Mar. 10) from 
the NOT (0.02A px -1 ). The profiles are H /3 (dot-dash), H7 (dash), and H<5 (solid). The 
maximum shock velocity was >3000 km s -1 . The stellar velocity is -54 km s -1 , note 
the wind absorption at an expansion velocity of about 10 km s -1 that becomes optically 
thin as the outburst progresses. The broad wings are due to a shock from the ejecta and 
subsequent outflow in the Mira (see discussion and Shore et al. 201 la, b). 

otics in outburst. The subsequent passage of the shock and the onset of wind 
recombination renders the profiles more symmetric with time, much line the NOT 
observation of H 7 and H<5 in AG Dra (Fig. 5b). 
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